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Abstract. Brueckner theory is used to investigate the properties of hyperons in nuclear matter. The 
hyperon-nucleon interaction is taken from chiral effective field theory at next-to-leading order with SU(3) 
symmetric low-energy constants. Furthermore, the underlying nucleon-nucleon interaction is also derived 
within chiral effective field theory. We present the single-particle potentials of A and 17 hyperons in sym¬ 
metric and asymmetric nuclear matter computed with the continuous choice for intermediate spectra. The 
results are in good agreement with the empirical information. In particular, our calculation gives a repulsive 
U-nuclear potential and a weak A-nuclear spin-orbit force. 
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1 Introduction 

The interaction between hyperons (Y) and nucleons (JV) 
is not only of interest by itself but it constitutes also the 
input for microscopic calculations of few- and many-body 
systems involving strangeness, such as hypernuclei or ex¬ 
otic neutron star matter. Indeed, with regard to the latter, 
the observation of two-solar-mass neutron stars m pro¬ 
vides strong restrictions for the appearance of hyperons in 
neutron star matter or, more generally speaking, on the 
in-medium properties of the hyperon-nucleon (YN) inter¬ 
action [31I31[5]. In particular, a sufficiently stiff equation-of- 
state (EoS) is required which does not leave much room 
for the presence of hyperons in the dense cores of neu¬ 
tron stars. A naive introduction of A-hyperons as an ad¬ 
ditional baryonic degree of freedom would soften the EoS 
such that it fails to support two-solar-mass neutron stars 
[^. There emerges a quest for repulsive A-nuclear forces 
at high baryon densities. Purely phenomenological mech¬ 
anisms have so far been invoked, e.g. through ad-hoc vec¬ 
tor meson exchange urn, multi-Pomeron exchange 0 or 
a suitably adjusted repulsive ANN three-body interaction 
[3]. Thus, there is an obvious need for a more systematic 
investigation of this issue. 

Various phenomenological approaches based on meson- 
exchange models [TUl[Tn[T^13Lll4j or quark models [TSlfTB] 
have been used to construct YN interactions in the past. 
However, given the poor experimental data base, these in¬ 
teractions differ considerably from each other. At the basic 
level the baryon-baryon interaction arises from the funda¬ 


mental theory of the strong interaction, quantum chromo¬ 
dynamics (QCD). Lattice QCD is approaching this prob¬ 
lem via numerical simulations [iiiiiHiiin]- Alternatively, 
chiral effective field theory (chiral EFT) is exploiting the 
symmetries of QCD together with the appropriate low- 
energy degrees of freedom to construct the baryon-baryon 
interactions. The YN interaction has been investigated at 
leading order (LO) [20] as well as next-to-leading order 
(NLO) pni^ within SU(3) chiral EFT. The YY inter¬ 
action between all members of the baryon octet has also 
been studied [2^l241l251l26j . The SU(2) chiral EFT frame¬ 
work, very successful in describing the low-energy nuclear 
forces to high precision P71I2H] . has thus been extended 
to the strangeness sector. Essential features of chiral EFT 
are that the results can be improved systematically by go¬ 
ing to higher order in the power counting scheme, and that 
two- and three-baryon forces can be calculated in a consis¬ 
tent way. These properties make chiral EFT very suitable 
for describing baryonic forces. The next-to-leading order 
YN potentials as derived from SU (3) chiral EFT include 
one- and two-pseudoscalar-meson (tt, K, rf) exchange dia¬ 
grams and four-baryon contact terms with SU(3) symmet¬ 
ric low-energy constants (LECs). Within this approach an 
excellent description of the experimental scattering data, 
comparable to the best phenomenological models, has been 
achieved m- 

Numerous advanced techniques have been developed 
to investigate systems without and with strangeness us¬ 
ing such microscopic interactions. For instance in the few- 
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body sector Faddeev-Yakubovsky theory [291130] can pro¬ 
vide very accurate results for systems with three or four 
particles. Many-body approaches such as Quantum Monte 
Carlo calculations |311l321l3^ . or nuclear lattice computa¬ 
tions IMlISaiSSj work well in the nucleonic sector and can 
be extended to the strangeness sector. Also many-body 
perturbation theory with chiral low-momentum interac¬ 
tions yields a good description of nuclear mat¬ 

ter. Concerning A and E hyperons in nuclear matter, spe¬ 
cific long-range processes related to two-pion exchange be¬ 
tween hyperons and nucleons in the nuclear medium have 
been studied in refs. gniiiT]. 

Conventional Brueckner theory [^I^l4^ at first order 
in the hole-line expansion, the so-called Bruecker-Hartree- 
Fock approximation, has been widely applied to calcula¬ 
tions of hypernuclear matter [ni[I3l43|46] with the use 
of phenomenological two-body potentials. The same ap¬ 
proach is also employed in the context of neutron star mat¬ 
ter |T7ll^l4^ . Recently, corresponding calculations of the 
properties of hyperons in nuclear matter have been also 
performed with a chiral YN interaction m- It has been 
demonstrated that the resulting in-medium properties of 
hyperons in symmetric nuclear matter around saturation 
density are consistent with the empirical information from 
hypernuclear phenomenology. In particular, a repulsive E- 
nuclear mean field and a very weak A-nuclear spin-orbit 
force have been found. 

In the present work we extend the hypernuclear mat¬ 
ter calculation of ref. [SD] in several ways: (i) We employ 
an underlying nucleon-nucleon interaction that originates 
likewise from chiral EFT, (ii) we implement the continuous 
choice for intermediate-state spectra in Brueckner ladders, 
and (iii) we investigate isospin-asymmetric nuclear matter 
with variable neutron and proton densities. Moreover, we 
compare our results to those obtained with the gap choice 
and with phenomenological YN potentials. Our calcula¬ 
tion of hyperons in infinite symmetric and asymmetric nu¬ 
clear matter should be considered as starting point for fu¬ 
ture studies of hypernuclei and exotic neutron star matter 
within SU(3) chiral EFT. 

Genuine three-baryon forces are disregarded in this 
work since the focus is on the hyperon-nucleon two-body 
interaction with special reference to its momentum de¬ 
pendence. At NLO in the chiral power counting no three- 
baryon forces arise. Chiral three-nucleon forces are, how¬ 
ever, very important in order to get saturation of nuclear 
matter from chiral low-momentum two-body interactions 
treated in many-body perturbation theory [3H]. The role 
of hyperon-nucleon-nucleon three-body forces for hyper¬ 
nuclei and exotic neutron star matter is a topic still under 
discussion, cf. refs. [41IM] . This will be subject of future 
studies within the applied chiral effective field theory. 

The paper is organized as follows. In sect, j^the con¬ 
struction of the baryon-baryon interaction from chiral ef¬ 
fective field theory is reviewed and the employed non- 
relativistic Brueckner-Hartree-Fock method is described 
in some detail. The results for hyperon single-particle po¬ 
tentials in isospin-symmetric and asymmetric nuclear mat- 



Fig. 1. Leading and next-to-leading order diagrams con¬ 
tributing to the baryon-baryon interaction potential. Solid 
and dashed lines denote octet baryons {N, A, E) and mesons 
{'K,K,ri), respectively. 


ter, as well as pure neutron matter, are presented in sect.j^ 
Section gives a summary of our findings. 


2 Formalism 

2.1 Baryon-baryon interaction 

For the description of the hyperon-nucleon interaction we 
use SU(3) chiral effective field theory up to next-to-leading 
order with the Weinberg power counting applied to the 
potential, as reported in detail in mm- At LO one- 
meson exchange diagrams and non-derivative four-baryon 
contact terms contribute to the potential. At NLO addi¬ 
tional contact terms and two-meson exchange diagrams at 
the one-loop level arise, cf. fig. The contact terms rep¬ 
resent the unresolved short-distance dynamics, and the 
corresponding low-energy constants (LECs) are fitted to 
low-energy YN scattering data and the hypertriton (J 4 H) 
binding energy. Given the sparse experimental information 
(in total 35 cross-section data points and a capture ratio) 
SU(3) flavor symmetric LEGs have been used, while SU(3) 
symmetry breaking is incorporated through the physical 
masses of the exchanged pseudoscalar mesons (tt, AT, ry). 
Additional constraints on the P-wave potentials come from 
the nucleonic sector via SU(3) symmetry. 

In order to obtain the reaction matrix (or T-matrix), 
the (two-particle irreducible) potentials are inserted into a 
regularized Lippmann-Schwinger equation, which involves 
coupled partial waves as well as coupled two-baryon chan¬ 
nels. The coupled-channel Lippmann-Schwinger equation 
in the particle basis reads after partial-wave decomposi¬ 
tion 


r°° dfcP 


E 


(27r) 


3 V^„^'\k",k) 


( 1 ) 


where J denotes the conserved total angular momentum, 
ly labels the two-particle channels {Ap, E^n, E^p ,...), and 
p represents the partial waves (^S'o, ■ ■ ■)■ Eurthermore 
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Fig. 2. Nucleon-nucleon phase shifts for the ^So and partial wave. The band represents the variation of our NLO results 
with the cutoff, see text. The black line is obtained with a cutoff A = 600 MeV. The circles denote the results from the GWU 
single-energy np partial wave analysis |53| . 


Table 1. Contact terms and threshold parameters for the ^So 
and NN partial waves for various cutoffs. The val¬ 

ues of the C”s are in 10"* GeV“^, the ones of the C”s, in 10^ 
GeV“^; the values of A in MeV. The scattering length a and 
the effective range r are in fm, the deuteron binding energy 
Ed in MeV. The empirical values are a-iso = —23.739 fm, 
ri_g^ = 2.68 fm, = 5.420 fm, = 1.753 fm |52| and 
Ed = -2.224575(9) MeV. 


A 

500 

550 

600 

650 

A27 

/^27 

*-"331 

^10* 

'-"3Si 

-0.1539 

2.313 

-0.2100 

0.2977 

-0.2767 

-0.1017 

2.326 

-0.1493 

0.3139 

-0.2896 

-0.0153 

2.326 

0.0166 

0.5109 

-0.2422 

0.1301 

2.328 

0.2059 

0.4899 

-0.2234 

“iSo 

-23.8 

-23.8 

-23.8 

-23.7 


2.81 

2.75 

2.68 

2.62 

“3 Si 

5.42 

5.43 

5.42 

5.43 

»'3Sl 

1.81 

1.76 

1.72 

1.67 

Ed 

-2.257 

-2.213 

-2.193 

-2.145 


^/s is the center-of-mass energy and stands for the per¬ 
tinent reduced baryon mass. The on-shell momentum ki, is 
determined by ^/s = ^ -|- k'^ + + k^. Note 

that high-energy components of the potential are cut off 
by a regulator function of the form fniA) = exp[—(fc'"* + 
fc^)/7l^], as done likewise in the nucleonic sector (54]. In 
the following the cutoff A is varied in the range of (500 
... 650) MeV for NLO and (550 ... 700) MeV for LO, i.e. 
in a range similar to that adopted for the nucleon-nucleon 
(NN) interaction in ref. [S3]. The resulting bands give an 
indication for the cutoff dependence and thus provide a 


lower bound on the theoretical uncertainty. Recently, im¬ 
proved schemes to estimate the theoretical error were pro¬ 
posed and applied to the NN interaction [ssiismisTj . How¬ 
ever, one has to go to higher orders in the chiral expansion 
if one wants to employ these schemes in a meaningful way. 
Therefore, we refrain from performing such an uncertainty 
quantification for the NLO interactions employed in the 
present study. 

A simultaneous description of the NN and YN in¬ 
teractions with SU(3) symmetric LECs is not possible at 
NLO |2T] (due to the strong correlation between the 
partial wave in the NN (/ = 1) and EN (7 = 3/2) chan¬ 
nels imposed by SU(3) symmetry). Therefore, we use dif¬ 
ferent sets of LECs in the NN and YN sectors. For the 
Y N interaction the set given in refs. |2T] and [SD] is usecQ 
where in the latter work the contact term that gives rise 
to an antisymmetric spin-orbit force in the YN interac¬ 
tion, responsible for spin singlet-triplet transitions, is al¬ 
ready fitted to the weak A-nuclear spin-orbit interaction 
msM- The NN interaction is based on the same dia¬ 
grams and contact terms as given in m, but with dif¬ 
ferent LECs, compiled in table [l] Furthermore, table [l] 
provides the scattering length, the effective range and the 
deuteron binding energy, calculated from this chiral NLO 
NN interaction. In fig. the nucleon-nucleon phase shifts 
in the and ^Si partial waves are shown. Note that the 
results are comparable to those (at NLO) in ref. [SD] where 
SU(2) chiral EFT was used (cf. fig. 4 in that reference). 
At low energies they are in agreement with the empirical 
data. However, at higher energies the results of the NLO 
interaction become too repulsive. 


^ Note that in order to be consistent with the definitions in 
eq. (18) of ref. | 50| the constants (7®“®“ in Table 1 have to be 
multiplied with a factor 2. 
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2.2 First order conventional Brueckner theory 

In order to investigate the properties of hyperons in nu¬ 
clear matter we employ the conventional Brueckner theory 
at first order in the hole-line expansion the so- 

called Brueckner-Hartree-Fock approximation. We focus 
on the single-particle potentials, i.e. the mean fields of hy¬ 
perons in nuclear matter. Our calculations are done in the 
particle basis and in the following we summarize the rele¬ 
vant formalism. For a more detailed introduction we refer 
the reader to ref. m and also to refs. [Tniiiiiniiig. 

A central object of Brueckner theory is the G-matrix 
representing the effective interaction between two parti¬ 
cles in the medium. The Brueckner reaction matrix or 
G-matrix is determined by solving the Bethe-Goldstone 
equation 

with e{u}) the energy denominator depending on the start¬ 
ing energy w. The Pauli blocking operator Q ensures that 
particles in intermediate states are outside their own Fermi 
sea. The same potential V as in the Lippmann-Schwinger 
equation Q for free scattering is used. Additional medium 
modifications of the two-meson exchange potential are 
not included. These represent density dependent two-body 
forces arising from three-body interactions. The medium 
effects come, therefore, solely from the Pauli operator Q in 
the Bethe-Goldstone equation and the density-dependent 
single-particle potential in the energy denominator e(a;). 

After angle-averaging, the Bethe-Goldstone equation 
decomposes into partial waves and reads for conserved val¬ 
ues of the total angular momentum J, total momentum 
K and starting energy w: 


Gffyik", k'; K,UJ) = k') 


^dfc k^ 




(27 


K — ki + k2 , k = 




Q,y{K,k) = - J dcos0 G 




( 2 ) 


^21 kp 


(5) 


when averaged over the angle 9 between K and fc, with the 
shorthand notation [a|6|c] = max(a, min(&, c)) introduced 
in [1^, and the arguments 


2:1 = 


22 = 



+k^-{k^pr 


+k^-{k^p^f 


( 6 ) 


The angle-averaged energy denominator is given by 

-ReC7B,(fci)-ReGB,(fc2), (7) 

with total and reduced masses, = Mi + M 2 and /i^ = 
M 1 M 2 /{Ml + M 2 ). The angle-average is done for the ar¬ 
guments of the single particle potentials t/s- of the inter¬ 
mediate baryons: 


ki = 


((i+Lf 


rAT^ + k ^+2 


i-H?i 


-Kk cos 6 


tional cosine 
f^j^dcosO cos 9 Q(K, k) 
J^^dcos 9 Q(K,k) 


1/1 


1/2 


( 8 ) 


with the mean directional cosine 

_ f' 

cos9 = 


= 2([-i|^2|i]-[-i|zi|i]) 


(9) 


where Q{K, k) is the exact Pauli blocking operator. If two 
nucleons are involved, the previous expression for cos 9 
would vanish (in symmetric nuclear matter), because of 
the almost equal masses, zi Ki Z 2 = zq. Then the alterna¬ 
tive angular average 


^ Gt'-'{k,k'-,K,ui). (3) 

e^(A, k; w) -I- le 

As in eq. Q the symbol p stands for the partial waves, 
p = {SL). The (coupled) two-particle channels are v = 
{B 1 B 2 ), with the baryons from the set: Bi G {n,p, A, 17+, 
i7°, 17“}. In the initial state the baryon B 2 is within its 
own Fermi sea. We introduce the total and relative mo¬ 
menta of two baryons Bi and B 2 by 

^ 12^1 — ^2 


COS0 = 


d cos 9 cos^ 9 Q(K, k) 1 

^ = "[o|^q|i] ( 10 ) 

\ f^^dcos 9 Q(K,k) 


is often used. 

It is common practice to introduce a further simplifi¬ 
cation. The squared momenta = K^(ki,k) and = 
fc|(A;i,A;) entering the Bethe-Goldstone equation are re¬ 
placed by their angle averages: 


K\ki,k) = 


J\kPCk(^)dcosd K^{ki,k,cosd) 


In eq. <© we have applied the standard approxima¬ 
tion replacing Q/ehy the ratio of its angle-averages Qje. 
The Pauli blocking operator involving the Fermi momenta 
( 12 ) 

kp' ' oi the two baryon species is given by 


= (1 + ^ 12 )^ [kf + k'^ — kik{l + [—Ijxoll])] , 


kl{ki,k) = 


ii2 


1 + ^12 


K^iki,k) + {l+^l2)k^ -^I2kf 


( 11 ) 


where 9 is the angle between ki and fc, and Xq means: 

^( 2^1 + (1 + Ci2)^ft2 - {kP)'^ 


Xo = 


2 ^ 12(1 + ^i 2 )kik 


( 12 ) 
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Note again, that the baryon B 2 in the initial state is within 
its Fermi sea. 

Finally, the single-particle potential of a baryon Bi 
due to the Fermi sea of the species B 2 is calculated in the 
Brueckner-Hartree-Fock approximation as follows 






x^(2J+l) 

J,P 


r-femax ^2 
fc ■ (27r)3 

^min ' ' 


X W{ki,k) 


(13) 


The full single-particle potential of a baryon Bi {Bi = 
n,p,A, is then given by the sum of the contributions 

from all baryons B 2 {B 2 = n,p) in the nuclear Fermi sea. 
The weight function W{ki, k) has the form 

VF(fci, fc) = ]- [ da = J(1 - [-l|xo|l]). (14) 

The lower and upper integration boundaries of the rel¬ 
ative momentum, kj^in and fcmax, are determined from 
W{ki,k) = 0, which leads to 


/cmin = niax 0, 


—kp'^ + ^i2ki 
1 + ^12 


a( 2 ) 


k — 

^ma.x — 


+ Cl 2^1 


1 + C 12 


(15) 


The G-matrix elements in eq. (13) are calculated at the 
on-shell starting energy 

Wo.s. = Eb^ (fci) -I- Eb 2 ih) 

Ebi (ki) = Mi + 2M- {ki). (16) 

This makes the determination of the single-particle po¬ 
tentials dependent on the si ngle -particle potential itself, 
and, therefore, eqs. (§ and ( [ 1 ^ have to be solved self- 
consistently. 

In the so-called gap choice the single-particle poten- 

di) 


for ki > kp\ Therefore, only the free particle energies 


tial is given by eq. (13) for ki < kp' and set to zero 


of the intermediate states appear in the energy denom¬ 
inator of the Bethe-Goldstone equation since the Pauli 
blocking operator is zero for momenta below the Fermi 
momentum. In the so-called continuous choice eq. (13) is 


used for the whole momentum range, hence the single¬ 
particle potentials enter also in the energy denominator. 
The result of the Brueckner theory to all orders in the 
hole-line expansion should be independent of the choice 
for the intermediate-state spectrum, however, at leading 
order the results differ somewhat. In the nucleonic sector 
at first order the continuous choice lies actually closer to 
the second order result of Brueckner theory, according to 
ref. [61j . Moreover, the continuous choice for intermediate 
spectra allows for a reliable determination of the single¬ 
particle potentials including their imaginary parts |45] . 
Unless stated differently, we use in this work the continu¬ 
ous choice. 


3 Results 


In this section we present our results for the in-medium 
properties of hyperons, based on a. YN interaction de¬ 
rived from chiral EFT [^. Additionally, for the ease of 
comparison, the G-matrix results obtained with two phe¬ 
nomenological YN potentials, namely of the Jiilich ’04 
m and the Nijmegen NSC97f [T^ meson-exchange mod¬ 
els, are given. Note that, like the EFT potentials, those 
YN interactions produce a bound hypertriton |30j . 

As mentioned before, the EFT NN and YN inter¬ 
actions involve different sets of low-energy constants. For 
calculations with the LO and NLO hyperon-nucleon inter¬ 
action we employ as the underlying nucleon-nucleon inter¬ 
action the NLO version with the same cutoff. In the case 
of the phenomenological YN interactions (Jiilich ’04 and 
Nijmegen NSG97f) we use for the purpose of comparison 
the NLO chiral NN potential with a cutoff of 600 MeV. 
In all calculations sums over partial waves up to J=5 are 
performed. 

First, we review the results for the nucleon single¬ 
particle potential derived from chiral effective field theory 
at NLO, as this is an input for calculations of hyperons in 
nuclear matter. Figure shows Uiq{k) for symmetric nu¬ 
clear matter at the Fermi momenta kp = 1.35 fm~^ and 
kp = 1.0 fm“^, corresponding to densities p = 0.166 fm“^ 
and p = 0.068 fm“^, as determined from the chiral NN 
potential and from the Nijmegen 93 model [52]. According 
to the Hugenholtz-van-Hove theorem the value at k = kp 
has to be Upfikp) « —53 MeV at saturation density. The 
results of our calculation with the EFT interaction are 
consistent with this constraint. Furthermore, fig. shows 
the total binding energy per particle 


E _ e 
A p ■ 


(17) 


The particle density p is given by a sum over the baryonic 
species that occupy Fermi seas. 


^ = E 

B 



(18) 


and the energy density e can be calculated from the single¬ 
particle potential as 


= E 


( kf^^ 

\ IOtt^Mb 


1 

27r2 


t(B) 


dkk^ ReUB{k) ) . (19) 


As it is typical for non-relativistic G-matrix calculations 
with realistic interactions the empirical saturation point of 
isospin-symmetric nuclear matter is not reproduced with¬ 
out the inclusion of three-nucleon forces [551151] . Note that 
the employed nucleon-nucleon interaction at NLO in chiral 
EFT becomes too repulsive for higher energies (cf. hg.j^. 
This feature appears to be reflected in the curve for the 
binding energy per nucleon which saturates at lower den¬ 
sities than usually found in calculations using (chiral and 
other) nucleon-nucleon potentials |641I651[661I67LI68| . But, 
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Fig. 3. Nucleon single-particle potential in symmetric nuclear matter for kp = 1.35 fm“^ and kp = 1.00 fm“^ (left) and energy 
per particle of nuclear matter with different proton fractions Pp/p (right). The bands represent the variation of our results with 
the cutoff, see text. The diamond symbolizes the empirical saturation point of symmetric nuclear matter. 


Table 2. A single-particle potential f7yi(fc = 0) in symmetric nuclear matter at saturation density, kp = 1.35 fm ^. Values are 
given in MeV and decomposed into partial wave contributions. 


UA{k = 0) 



"Fo 

^Pi 



Total 

NLO (500) cont 

-15.4 

-15.7 

1.0 

1.8 

1.5 

-1.3 

-28.3 

NLO (550) cont 

-13.9 

-12.7 

0.9 

1.6 

1.5 

-1.2 

-24.2 

NLO (600) cont 

-12.9 

-13.5 

0.8 

1.3 

1.4 

-1.2 

-24.4 

NLO (650) cont 

-12.4 

-16.3 

0.7 

1.2 

1.3 

-1.2 

-27.0 

LO (600) gap 

-12.1 

-25.9 

-1.7 

1.5 

1.7 

-0.4 

-37.2 

LO (600) cont 

-13.2 

-28.0 

-1.9 

1.5 

1.7 

-0.4 

-40.7 

NLO (600) gap 

-13.1 

-13.9 

0.9 

1.3 

1.4 

-1.2 

-24.8 

NSC97f gap 

-14.7 

-24.1 

0.4 

2.4 

4.1 

-0.8 

-34.1 

NSC97f cont 

-14.5 

-25.2 

0.4 

2.3 

3.9 

-0.9 

-35.5 

Jiilich ’04 gap 

-10.5 

-36.5 

-0.7 

-0.6 

0.5 

-3.2 

-51.7 

Jiilich ’04 cont 

-11.2 

-38.0 

-0.7 

-0.7 

0.5 

-3.3 

-54.2 


as expected, our results still lie within the well-known Co¬ 
ester band [69]. In this context we want to stress that we 
show the nucleonic results only for illustrative purposes. 
Considering the recent arrival oi NN interactions at fifth 
order in chiral EFT izniin] the NLO potential employed 
here is obviously not state-of-the-art. However, for consis¬ 
tency reasons we prefer to use NN and YN interactions at 
the same order of the chiral expansion. In any case, as we 
will see below, the properties of hyperons in nuclear mat¬ 
ter do not depend strongly on the nucleon single-particle 
potential U]\f{k), and therefore the NN interaction up to 
NLO is certainly sufficient for our purposes. 

Now we turn to the properties of hyperons in sym¬ 
metric nuclear matter as they follow from SU(3) chiral 
effective field theory. In table [Rvalues for the depth of the 
A single-particle potential C/yi(/c = 0) at saturation den¬ 


sity are given. The results of the LO and NLO results are 
consistent with the empirical value of about —28 MeV as 
deduced from binding energies of A hypernuclei 
The results for the gap choice of intermediate spectra are 
similar to ref. |50j , where a phenomenological parametriza- 
tion of the nucleon single-particle potential Lbv(^) has 
been used. This suggests that the results for the A single¬ 
particle potential do not depend strongly on those of the 
nucleon in nuclear matter. (Actually, using this phenomeno¬ 
logical parametrization, we can reproduce the results of 
ref. (SD) which served as a test for the new code devel¬ 
oped for the present investigation.) The differences for 
UA{k = 0) between the gap choice and the continuous 
choice are a few MeV, comparable to what has been found, 
e.g., in ref. [12] for the Nijmegen NSC97 potentials. Ob¬ 
viously, the two phenomenological models (Jiilich ’04, Ni- 
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Fig. 4. Momentum dependence of the real and imaginary parts of the single-particle potential of a 4 hyperon in isospin- 
symmetric nuclear matter at saturation density. The bands represent the variation of our results with the cutoff, see text. 

Table 3. S single-particle potential Us{k — 0) in symmetric nuclear matter at saturation density, kp = 1.35 fm~^. Values are 
given in MeV and decomposed into partial wave contributions. 


Usik = 0) 


^Si+^Di 

"Po 

^Pi 

"Pi 

"P2+"P2 

Total 

NLO (500) cont 

-4.6 

13.5 

1.5 

0.6 

0.4 

0.3 

11.6 

NLO (550) cont 

-4.2 

17.3 

1.3 

0.4 

0.2 

-0.0 

14.9 

NLO (600) cont 

-4.7 

15.4 

1.2 

0.2 

0.0 

-0.4 

11.5 

NLO (650) cont 

-4.9 

11.9 

1.0 

0.0 

-0.1 

-0.8 

7.0 

LO (600) gap 

-1.8 

25.3 

-1.9 

-0.2 

-1.4 

-1.1 

18.7 

LO (600) cont 

-2.2 

22.1 

-1.9 

-0.2 

-1.2 

-1.0 

15.5 

NLO (600) gap 

-5.6 

15.4 

1.1 

0.1 

-0.2 

-0.6 

9.9 

NSC97f gap 

1.9 

-17.3 

0.4 

-2.1 

1.1 

-2.4 

-19.1 

NSC97f cont 

-0.5 

-17.2 

0.4 

-2.1 

1.0 

-2.6 

-22.0 

Jiilich ’04 gap 

-8.4 

-4.0 

0.4 

-2.0 

-1.8 

-3.7 

-20.1 

Jiilich ’04 cont 

-7.9 

-5.6 

0.4 

-2.0 

-2.5 

-3.8 

-21.9 


jmegen NSC97f) predict more attractive values of Ua{0) = 
(—35 .. .—50) MeV, where the main difference is due to the 
contribution in the ^Si partial wave. As already discussed 
in ref. [SD], we believe that this due to the fact that the 
Ap transition is significantly larger in the NLO 

chiral EFT interaction as compared to the one of the LO 
interaction and of the phenomenological models whereas 
the diagonal and transitions are ac¬ 

cordingly smaller. 

The momentum dependence of the real and imagi¬ 
nary parts of the A single-particle potential is presented in 
fig-11 A marked difference between LO and NLO is, that 
the A single-particle potential at NLO turns to repulsion 
at fairly low momenta around fc « 2 fm~^. A similar be¬ 
havior is also found for the NSC97f potential. The cutoff 
dependence at LO seems to be accidentally weak. 


Corresponding results for E hyperons in isospin-sym- 
metric nuclear matter at saturation density are given in 
tableland are also graphically displayed in fig.[^ The pre¬ 
sented results are for the neutral E^ hyperon. The small 
difference to the results for charged E^ and E~ hyper¬ 
ons comes solely from the mass difference of the three 
E particles, and is of the order of (0.5... 1) MeV, where 
the difference between and E^ is larger than the one 
between and E~. According to analyses of data on 
{tt~ , K^) spectra related to E~ formation in heavy nuclei 
the V-nuclear potential is moderately repulsive in sym¬ 
metric nuclear matter, see the review |74] . This feature is 
well reproduced in our calculation for NLO and even for 
LO. Indeed, in the course of constructing the NLO inter¬ 
action it turned out that the available YN scattering data 
could be fitted equally well with an attractive or a repul¬ 
sive interaction in the ^Si partial wave of the I = 3/2 EN 
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Fig. 5. Momentum dependence of the real and imaginary parts of the single-particle potential of a If hyperon in isospin- 
symmetric nuclear matter at saturation density. The bands represent the variation of our results with the cutoff, see text. 


channel [H], which is the partial wave that provides the 
dominant contribution to the S single-particle potential, 
cf. tablej^and also table 4 in [50] . For the reasons discussed 
above, the repulsive solution was adopted. Note that mod¬ 
els derived within the meson-exchange framework often 
fail to produce a repulsive if-nuclear potential and the 
two phenomenological YN potentials considered here are 
exemplary for this deficiency. As visible in fig. the S 
potential stays repulsive for higher momenta. The imagi¬ 
nary part of the A-nuclear potential at saturation density 
is in good agreement with the empirical value of —16 MeV 
as extracted from A”-atom data m- The imaginary po¬ 
tential is mainly induced by the EN to AN conversion 
in nuclear matter. Evidently, the bands representing the 
cutoff dependence of the chiral potentials, become smaller 
when going to higher order in the chiral expansion. This 
feature has been also observed for the YN scattering ob¬ 
servables m- 

In the following we provide a more detailed view on 
the dependence of our in-medium results on the densities 
of protons and neutrons. The corresponding predictions 
are shown only for the chiral EFT interaction at NLO 
with a fixed cutoff, namely A = 600 MeV, for reasons of 
clearer presentation. However, one should keep in mind 
that these results are likewise subject to variations with 
the cutoff and, specifically, in case of a weak sensitivity 
to the density the latter effect might be actually smaller 
than the cutoff dependence. 

In figs. 1^ and the dependence of the single-particle 
potential C/y(fc) for A and E on the density in symmet¬ 
ric nuclear matter is shown. The chosen nucleon Fermi 
momenta kp = (1.00,1.35,1.53) fm“^ correspond to den¬ 
sities of about p = (0.4,1.0, 1.5)/9o with po = 0.16 fm“^. 


The momentum dependence of the potentials is similar 
for different densities, but their magnitude varies strongly. 
Especially for the E hyperon the single-particle potential 
can even become attractive at low densities. 

Pure neutron matter is another interesting environ¬ 
ment for the in-medium behavior of hyperons. Therefore, 
we display in fig. also the density dependence of f7y(fc) 
for A and E hyperons in pure neutron matter. The Fermi 
momenta of the neutrons kp = (1.26,1.70,1.92) fm“^ cor¬ 
respond the same densities as selected in figs. and 
for isospin-symmetric nuclear matter. Due to the maxi¬ 
mal asymmetry between protons and neutrons, the single¬ 
particle potentials for the three {E~^, E^, E~) hyperons 
are rather different. The A single-particle potential UA{k) 
is slightly more attractive than the one in isospin-sym- 
metric nuclear matter at the same density, cf. fig. thus 
indicating only a weak dependence on the composition of 
nuclear matter. The small surplus of attraction can be 
understood from the reduction of Pauli blocking effects in 
pure neutron matter. 

In order to get a more detailed insight into the inter¬ 
action of hyperons with heavy nuclei, we consider also the 
strength of the A-nuclear spin-orbit coupling. It is experi¬ 
mentally well established cniiiTi that the A-nucleus spin- 
orbit force is very small. In the following we will present 
results for the so-called Scheerbaum factor It quanti- 
hes the strength of the A- or V-nuclear spin-orbit poten¬ 
tial, which takes the form [75] 


U^£{r) 



1 dp(r) 
r dr 


I ■ cr , 


( 20 ) 


where p{r) is the nucleon density distribution, I the (single¬ 
particle) orbital angular momentum operator and cr the 
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Fig. 6. Momentum dependence of the real and imaginary parts of the A single-particle potential for different Fermi momenta 
kp = (1.00,1.35,1.53) fm“^ in symmetric nuclear matter, calculated in xEFT at NLO with a cutoff A = 600 MeV. 




Fig. 7. Momentum dependence of the real and imaginary parts of the E single-particle potential for different Fermi momenta 
kp = (1.00,1.35,1.53) fm“^ in symmetric nuclear matter, calculated in xEFT at NLO with a cutoff A = 600 MeV. 


hyperon spin operator. In isospin-symmetric nuclear mat¬ 
ter with Fermi momentum kp the Scheerbaum factor 
is obtained from the G-matrix elements via the relation m 

X ^ ^(2J-|-1)^12(1 + Cl2)^ 

B 2 —n,p J 


X 



dk 


IF(fci,A:)Re| 


(27r)3 

/ j I ^U+ 1 ,U+ 1 ,J 


11 

^CBiB2)(BiB2 ) 


-(J-l)G 
-^J[J + l)G 
-y/j{J+l)G 


(BiSalfSiBa) 

U.OJ.J 

(BiB 2 )(BiB 2 ) 

OJ.IJ.J 

(BiB 2 )(BiB 2 ) 


{k,k;K,LUo.s.) 

{k,k;K,LUo.s.) 

{k,k;K,LUo.s.) 

{k,k;K,uJo.s.) 

(/c, fc;Rr,Wo.s.)|, ( 21 ) 
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Fig. 8. Momentum dependence of the real part of the single-particle potentials for hyperons in pure neutron matter for different 
Fermi momenta kp = (1.26,1.70,1.92) fm“^, calculated in xEFT at NLO with a cutoff A = 600 MeV. 


with ki set to zero in the end. In table we present the 
Scheerbaum factor Sa for the A hyperon in symmetric nu¬ 
clear matter at saturation density. The values in table 
are in agreement with the earlier results of ref. |5()j . The 
difference between the gap and the continuous choice for 
intermediate spectra is small. This is expected because 
the Scheerbaum factor involves only contributions from 
P-waves and higher partial waves and it is known that 
these are much less sensitive to the treatment of the inter¬ 
mediate spectra than the 5'-waves [l2]. We use the same 
YN interaction as in ref. [SD] where the strength of the 
antisymmetric spin-orbit contact interaction, generating a 
spin singlet-triplet mixing (^Pi -o- ^Pi), has been tuned 
to achieve Sa ~ —3.7 MeVfm®, in accordance with es¬ 
timates for the empirical value that is expected to be in 
the range of around —4.6 to —3.0 MeV fm® [801IM] . In ta¬ 
ble [S] we summarize our results for the Scheerbaum factor 
Ss, which are close to the values reported in ref. [50] . 
As before, results are only given for The difference 
among the Scheerbaum factors for the three E hyperons 
due to their mass splitting is smaller than 0.5 MeVfm®. 
In contrast to the leading-order approximation, at NLO 
negative values of are always obtained, similar to the 
results found with the NSC97f and Jiilich ’04 models. 

Another important quantity to characterize in-medium 
properties is the effective baryon mass. The ratio between 
the effective and the free hyperon mass in nuclear matter 
is usually defined as 


Mb 


1 + 2Mb 


dReUB{k) 


-1 


k-O- 


( 22 ) 


We follow ref. [T5j where this ratio is computed as 


Mh 

Mb 



2Mb 

-^Re{UB{k)-UBm 


(23) 



P/Po 


Fig. 9. Density dependence of the effective mass of a A hy¬ 
peron in (a)symmetric nuclear matter, calculated in yEFT at 
NLO with a cutoff A — 600 MeV. 


with A: « 1 fm“^. In some works k is set to the Fermi 
momentum kp in symmetric nuclear matter. In fig. the 
density dependence of the effective mass M\ of a A hy¬ 
peron for different isospin asymmetries Pp = (0,0.25, 0.5)p 
is shown, where pp is the proton density and p the total 
density. The effective A mass in pure neutron matter is 
slightly higher than the one in symmetric nuclear matter, 
but the shape of the curves in fig. does not depend on 
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Table 4. Scheerbaum factor S/i(fc = 0) in symmetric nuclear matter at saturation density, kp = 1.35 fm ^. Values are given in 
MeV fm® and decomposed into partial wave contributions. 


SAik = 0) 




iPi 

®F2 


"Fs 

Total 

NLO (500) cont 

-5.6 

-0.6 

-4.4 

10.4 

-3.5 

0.4 

0.2 

-3.0 

NLO (550) cont 

-4.9 

-0.6 

-4.2 

9.2 

-3.2 

0.4 

0.2 

-3.1 

NLO (600) cont 

-4.4 

-0.6 

-4.0 

8.3 

-3.1 

0.4 

0.2 

-3.2 

NLO (650) cont 

-4.0 

-0.6 

-3.8 

7.3 

-3.1 

0.4 

0.2 

-3.6 

LO (600) gap 

9.4 

-0.2 

-4.9 

0.0 

-1.1 

0.4 

-0.1 

3.5 

LO (600) cont 

10.2 

-0.2 

-4.7 

0.0 

-1.2 

0.4 

-0.1 

4.5 

NLO (600) gap 

-4.7 

-0.6 

-4.1 

8.4 

-3.1 

0.4 

0.2 

-3.4 

NSC97f gap 

-2.2 

0.5 

-11.4 

2.0 

-2.6 

1.1 

-2.0 

-14.2 

NSC97f cont 

-2.0 

0.5 

-10.8 

1.9 

-2.8 

1.2 

-2.0 

-13.8 

Julich ’04 gap 

4.0 

0.4 

-1.4 

5.1 

-9.1 

0.6 

-1.0 

-1.3 

Julich ’04 cont 

4.1 

0.5 

-1.3 

5.1 

-9.3 

0.6 

-1.1 

-1.4 


Table 5. Scheerbaum factor Ss{k = 0) in symmetric nuclear matter at saturation density, kp = 1.35 fm Values are given 
in MeV fm® and decomposed into partial wave contributions. 


Sp{k = 0) 



®Fi ’ 

Fi ®Fi 

®F2 



Total 

NLO (500) cont 

-7.7 

0.6 

-1.0 

-10.6 

0.9 

0.1 

0.1 

-17.7 

NLO (550) cont 

-6.6 

0.3 

-0.6 

-9.7 

-0.0 

0.1 

0.1 

-16.5 

NLO (600) cont 

-5.9 

0.1 

-0.1 

-9.2 

-1.0 

0.1 

0.0 

-16.0 

NLO (650) cont 

-5.3 

0.1 

0.2 

-8.3 

-2.0 

0.1 

-0.0 

-15.3 

LO (600) gap 

9.5 

0.3 

3.2 

0.0 

-2.8 

0.1 

-0.4 

9.8 

LO (600) cont 

8.7 

0.3 

2.6 

0.0 

-2.4 

0.1 

-0.3 

8.9 

NLO (600) gap 

-5.6 

0.1 

0.4 

-10.1 

-1.3 

0.1 

-0.1 

-16.4 

NSC97f gap 

-2.0 

0.1 

-2.8 

-2.7 

-6.2 

0.3 

-1.5 

-14.9 

NSC97f cont 

-1.8 

0.0 

-2.5 

-2.6 

-7.0 

0.3 

-1.6 

-15.5 

Julich ’04 gap 

-2.2 

0.1 

4.5 

-8.2 

-9.6 

0.3 

-1.0 

-16.4 

Julich ’04 cont 

-2.2 

0.1 

6.0 

-9.8 

-9.9 

0.3 

-1.0 

-16.9 


the composition of nuclear matter. Since the momentum- 
dependence of the if-nuclear potentials as obtained from 
the G-matrix calculations is not close to a quadratic be¬ 
havior for momenta 0 < k < 3 fm~^, the effective mass 
does not serve as a significant quantity. 

Finally, we investigate in more detail the influence of 
the composition and density of nuclear matter on the 
single-particle potentials of hyperons. In fig.[^ the den¬ 
sity dependence of the depth of the nuclear mean-held of 
Aor S hyperons at rest is shown for isospin-symmetric nu¬ 
clear matter, asymmetric nuclear matter with pp = 0.25p 
and pure neutron matter. The single-particle potential of 
the A hyperon is almost independent of the composition of 
the nuclear medium, because of its isosinglet nature. Fur¬ 
thermore, it is attractive over the whole considered range 
of density 0.5 < pjpo < 1.5. The three E hyperons pos¬ 
sess (up to small differences from the mass splittings) the 
same single-particle potential in symmetric nuclear mat¬ 
ter. It is attractive for low densities, but turns into re¬ 
pulsion at p « 0.7po and stays repulsive for higher densi¬ 
ties. When introducing isospin asymmetry in the nuclear 
medium a splitting of the single-particle potentials occurs 
due to the strong isospin dependence of the EN interac¬ 
tion. The splittings among the E'^, E^ and E~ potentials 
as obtained in our microscopic calculation have a non¬ 
linear dependence on the isospin asymmetry which goes 


beyond the usual (linear) parametrization in terms of an 
isovector Lane potential |5^ . 

4 Summary and outlook 

We have investigated the properties of hyperons in nuclear 
matter at variable density and isospin asymmetry within 
the self-consistent Brueckner-Hartree-Fock approach us¬ 
ing the continuous choice for intermediate spectra. The 
employed microscopic potentials for the hyperon-nucleon 
interaction are derived within SU(3) chiral effective held 
theory at next-to-leading order. In addition the relevant 
nucleon-nucleon interaction is taken as well from chiral 
EFT. We have presented results for the complex-valued 
single-particle potentials of A and E hyperons and have 
studied the density dependence of the hyperon potential 
depths and the A effective mass. Our results reveal that 
the hyperon-nucleon interaction in chiral effective held 
theory is constructed in such way that it is not only con¬ 
sistent with YN scattering data, but also in agreement 
with empirical properties of the hyperon-nuclear poten¬ 
tials in symmetric nuclear matter. In particular, the po¬ 
tential depth 1 / 4 ( 0 ) « —28 MeV of the A hyperon is re¬ 
produced and the E single-particle potential comes out 
repulsive in accordance with recent measurements of E~ 
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Fig. 10. Density dependence of the hyperon single-particle potentials at fc = 0 with different compositions of the nuclear matter, 
calculated in ^EFT at NLO with a cutoff A = 600 MeV. 


formation spectra on heavy nuclei. Furthermore, the ex¬ 
perimentally indicated very small value of the yl-nuclear 
spin orbit coupling is reproduced. 

The presented calculations provide a basis for studies 
of hypernuclei. An extension of the calculations to higher 
densities of about 2-3po will allow for a first exploratory 
study of neutron star matter including hyperons from the 
viewpoint of chiral effective field theory. In particular, the 
inclusion of three-baryon forces, especially the ANN in¬ 
teraction, represents a very interesting forthcoming ex¬ 
tension of these calculations which can now be performed 
consistently within chiral effective field theory since the 
SU(3) classification of the three-baryon forces has been 
completed [831184] . 

This work is supported in part by the DFG and the NSFC 
through funds provided to the Sino-German CRC 110 “Sym¬ 
metries and the Emergence of Structure in QGD”. 
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